Applied Soil Ecology

Biochar amendments are increasingly applied to improve soil physical functioning and support carbon dioxide removal, but their effects on intrinsic soil thermal properties remain poorly characterised. We conducted the first global systematic meta-analysis of 19 independent studies, 231 control-biochar comparisons, and 529 property-specific effect sizes to test how biochar changes soil heat transfer and storage. Biochar reduced thermal conductivity by 17.6% (95% CI, -22.7 to -12.2), thermal diffusivity by 11.0% (-14.5 to -7.3), and volumetric heat capacity by 8.3% (-12.3 to -4.1). Gravimetric heat capacity showed no significant overall response (+3.3%; -7.6 to 15.4) but was supported by fewer studies. Negative responses were directionally consistent for thermal conductivity, diffusivity, and volumetric heat capacity. Moderator analyses showed that responses were most consistently associated with post-application bulk density and changes in bulk density, while application rate modulated response magnitude and soil texture constrained context dependence. Co-variation among thermal conductivity, thermal diffusivity, and volumetric heat capacity matched expected physical dependencies, indicating coordinated structural reorganisation rather than independent shifts in isolated parameters. These estimates describe intrinsic conductive and storage properties; field-scale soil temperature responses may also be modified by albedo, evaporation, vegetation, and surface energy balance. Improved integration of soil thermal measurements with moisture dynamics, structural changes, and carbon cycling is essential to accurately represent biochar effects in soil and land-surface models.

Streptomyces are abundant soil inhabitants with extensive secondary metabolism and antibiotic resistance traits. Yet, their ecological role in shaping soil antibiotic resistome dynamics remains understudied. Here, we investigated how two different bio-based fertilizers harbouring Streptomyces shaped soil resistome and mobilome by combining genome analysis of eight Streptomyces isolates to metagenomic profiling of soils before fertilization, within 48 hours after fertilizer application, and six weeks after. Streptomyces genomes showed linkages among antibiotic resistance genes, carbohydrate-active enzymes, and antibiotic-production-associated biosynthetic gene clusters, connecting resistance and biosynthesis to broader metabolic strategies. Relationships between carbon degradation and biosynthesis associated with specific enzyme families, indicating that carbon availability shapes secondary metabolism. We confirmed experimentally that antibacterial potential varied with carbon source, suggesting that microbial activity during manufacturing of the bio-based fertilizers may create localized selection pressures before fertilizers enter the soil. Fertilization with the studied materials induced modest but consistent shifts in resistome and mobilome without major changes in dominant taxa or overall bacterial abundances, indicating functional reorganization within soil communities. Diversity of antibiotic resistance genes and mobile genetic elements increased, whereas abundance changes were small. Mobile genetic element composition showed stronger responses that were associated with fertilizer inputs, Streptomyces abundance, and taxa linked to faecal and resistance sources. Together, our results show that bio-based fertilizers shape soil resistome primarily through ecological restructuring of resident soil communities, while carbon-dependent microbial activity within fertilizers may enrich resistance. These factors should be considered in manufacturing of bio-based fertilizer as well as in designing agricultural practices.

Secondary forest succession following agricultural abandonment is a dominant land-use transition across the tropics, yet whether soil microbial communities recover toward old-growth forest reference states remains poorly understood, particularly in West Africa. Here, we investigated the successional dynamics of bacterial and arbuscular mycorrhizal (AM) fungal communities along post-agricultural chronosequences spanning 1 to 43 years across six classified forests in Cote dIvoire, using Bayesian hierarchical models applied to amplicon sequencing data. Both guilds attained moderate to high alpha diversity within the first decade of succession; AM fungal diversity showed moderate evidence of age-related increase thereafter while bacterial diversity showed no directional trend. Pairwise turnover analyses revealed progressive internal convergence in AM fungal communities with plots farther apart in successional time becoming more compositionally similar, while bacterial communities showed only a weak and uncertain tendency in the same direction. Beta-dispersion modelling further indicated progressive within-forest homogenisation of AM fungal communities across abundance-weighted metrics, while bacterial assemblages showed no such stabilisation. Despite this internal convergence, compositional distances to old-growth reference plots remained persistently high for both guilds throughout the chronosequence, with no statistical evidence of recovery toward old-growth states across any dissimilarity metric or guild within the 40-year window. Indicator species analysis identified no robust stage-specific taxa after correction for multiple testing. These results indicate that microbial succession in post-agricultural West African forests is characterised by rapid early reorganisation followed by stabilisation into site-specific assemblages that remain persistently distinct from old-growth reference communities. This outcome challenges the direct application of classical vegetation successional theory to soil microbiomes and suggests that passive regeneration alone is unlikely to restore old-growth microbial communities within restoration-relevant timescales.

Soil microbial diversity is shaped by the spatial scale at which communities are sampled, yet standard sampling practices often homogenize samples, obscuring fine-scale spatial structure and diversity patterns. To better understand how sampling effort, spatial extent, and physical homogenization influence plot-level microbial richness estimates, we sampled 57 forest and grassland sites across three regions in Germany using a 14-core cross-transect design and performed 16S rRNA gene metabarcoding. We simulated sampling efforts and a range of spatial extents and compared diversity estimates to those from physically homogenized composite samples. Plot-level richness increased continuously with sampling effort and spatial extent, with no evidence of saturation. However, when sequencing depth was held constant, sampling completeness declined with increasing sampling effort, meaning that more diversity is not captured. Composite samples substantially underestimated plot-level richness and altered apparent diversity relationships between ecosystems; individual cores identified forests as richer than grasslands, whereas homogenized samples suggested the opposite relationship. These results demonstrate that sampling effort, spatial extent, and homogenization fundamentally shape soil microbial diversity estimates. Homogenized composite samples cannot substitute for individual cores when the goal is to reliably quantify plot-level richness or compare diversity across ecosystems.

The contribution of predation between bacteria to microbial community dynamics in agricultural fields has hardly been investigated. Here. dynamics of general prokaryotes (GEP) and of the predators Bdellovibrionales (Bd) and Bacteriovoracales (Bac) (Bdellovibrio-and-Like Organisms, BALOs) were studied in two agricultural fields differing in organic and mineral input regimes, for one year. Season, but not fertilization, affected absolute sizes of GEP and of BALO communities. 16S rRNA gene community analysis identified numerous novel Bd and Bac lineages, with none of the dominant BALOs related to characterized isolates. A few dominant BALO amplicon sequence variants (ASVs) persisted year-round, whereas others showed seasonal- or treatment specific responses. GEP, Bd, and Bac ASV a-diversity was mostly influenced by season, with some changes due to fertilization in Bd, and Bac communities. Seasonal changes, site, and fertilization regimes influenced {beta}-diversity of GEP, Bd and Bac communities and determined the structure of BALO-gram-negative bacteria interaction networks, signaling that niche segregation acts at the microbiome-BALO interface. Accordingly, we suggest that shifts in GEP community structure triggered by environmental changes and agricultural practices cascade to BALO predators, in turn affecting BALO-microbiome interactions. These dynamics may be harnessed to manipulate the soil microbiome to benefit sustainable environmental and agricultural outcomes.

Plant roots are broadly colonized by endophytic fungi with saprotrophic capabilities, but our understanding of whether they function in ways that are beneficial or detrimental to the host remains limited to model organisms. We hypothesized that endophytic fungi broadly affect plant access to soil nutrients, particularly organic forms that are typically not directly available to the plant. To address this, we paired 41 fungal endophytes with switchgrass (Panicum virgatum L.) and provided either inorganic or organic forms of nitrogen (N) and phosphorus (P). We evaluated how the fungi affected plant tissue N and P as well as plant growth. We also examined if these outcomes could be predicted from fungal phylogenetic relationships, in vitro traits of the fungi, or characteristics of the habitat from which fungi were isolated. There was substantial variation in both plant N (0.05-0.63%) and P (0.02-0.10%) acquisition that depended on the interaction of fungus and nutrient treatment. More fungi were beneficial for plant N than for P and shoot nutrients generally increased more than root nutrients from fungal associations. However, fungal effects on plant nutrients were not predicted by fungal traits, habitat traits, or fungal phylogenetic relationships. This unpredictability highlights a key challenge for incorporating endophytes into nutrient management strategies. Improving our ability to predict endophyte impacts on host nutrient acquisition will require identifying the mechanisms underlying observed beneficial effects and scaling up to realistic, diverse root microbial communities.

Animal COI (mitochondrial cytochrome oxidase I) metabarcoding of environmental DNA (eDNA) is increasingly used to assess biodiversity in complex substrates such as soil. However, due to read-length constraints of second-generation sequencing platforms, mini-barcodes have been used instead of the full barcode region. Long-read sequencing technologies now enable the recovery of full-length barcode sequences, and are more commonly applied for studying microbes, but their use for metabarcoding the full-length standard COI barcoding region in animals remains limited. In this study, we compared three COI amplicon sets -- 313 bp, 660 bp, and 1,256 bp -- amplified from soil eDNA samples and sequenced using Illumina and PacBio platforms to evaluate their overall concurrence, the effectiveness of identifying nuclear mitochondrial DNA segments (NUMTs) and chimeras, as well as their respective taxonomic resolution. The long-read datasets exhibited a higher identification rate of NUMTs and true chimeras, suggesting that longer sequences improve the detection of noise in COI metabarcoding data, thereby reducing the occurrence of spurious taxa. Taxonomy assignment confidence was similar between the 313 bp and 660 bp datasets, whereas extending the amplicon beyond the standard COI barcode region (1,256 bp) reduced confidence, likely because longer reads extend into regions poorly represented in barcode reference databases. Despite substantially lower sequencing depth in the 660 bp dataset, per-sample OTU richness did not differ significantly from that recovered with the Illumina 313 bp amplicon set. Similarly, the relationships between samples were strongly correlated across the detected OTU communities, indicating consistent ecological interpretations between short and long amplicons. We conclude that the standard ~658 bp COI barcode is an optimal marker for soil animal metabarcoding from eDNA, balancing target recovery, artifact detection, taxonomic assignment and ecological interpretability. As COI eDNA metabarcoding becomes increasingly used in biodiversity assessment and is increasingly adopted in large-scale monitoring initiatives, this study provides methodological guidance for improving the robustness of soil animal community biomonitoring.

Soils are recognized as reservoirs of antibiotic resistance genes (ARGs) with the potential to transfer to clinical pathogens, creating antimicrobial resistance (AMR) that poses a threat to human health. While large-scale AMR surveys have profiled how diverse biomes shape soil resistomes, less is known about the influence of specific soil properties. Here, we combined metagenomics and 16S rRNA amplicon sequencing with isolate-based approaches to investigate drivers of soil AMR across a tropical watershed from beach to mountaintop in Waimea Valley, Oahu, Hawai{square}i. We leveraged functional- and taxonomic-classification of resistances to unravel how soil properties interact with bacterial taxa to structure resistomes. Metagenomic- and isolate-resistomes showed remarkable consistency, including a general gradient of increasing AMR from ridge to beach. Resistome functional composition was significantly correlated with total bacterial community structure. The relationship between resistances and soil properties was primarily dictated by taxonomic composition of each resistance. Rifampin- and Vancomycin-ARGs associated with Actinomycetes negatively correlated with soil physical properties, while resistant genes and isolates from Gammaproteobacteria positively correlated with enzymatic activity metrics. These findings indicate that soil properties structure the resistome indirectly through taxonomic filtering of microbial hosts and challenge the notion that AMR is decoupled from phylogenetic relatedness.

Plants form holobionts by associating with diverse microbiota. Self-organization gives rise to emergent properties of the holobiont, such as increased resistance to pathogens. However, the local factors contributing to the self-organization are not well understood. We hypothesized that nematode communities and their associated microbiota govern the rhizobiome of the model plant holobiont tomato in terms of its suppression of root invasion by the parasite Meloidogyne hapla, and that the soil legacy influences the suppressive potential mediated by these biota. In pot experiments, a resistant tomato holobiont was favored by assembly in the presence of a nematode community conditioned by tomato plants, compared to oilseed rape or fallow soil. Nematode communities conditioned by tagetes could enhance resistance even better than tomato. Microbiota from crushed tomato-conditioned nematode communities increased resistance of the tomato holobiont, compared to microbiota of nematode communities conditioned by maize, or heat-inactivated microbiota. The 0.2 micrometre filtered microbiota from crushed nematodes had the same effect, suggesting a role of nematode-associated bacteriophages in holobiont assembly. The results indicate that soil nematodes and their associated microbiota play a role in the local organization and stabilization of plant holobionts. They can influence the resistance of plants that subsequently grow in the same soil. From an applied perspective, crop rotation schemes that alter nematode-microbiota communities could be harnessed to engineer crop holobionts.

Understanding the relationships between soil macrofauna and decomposition is crucial for predicting how land-use change impacts ecosystem function in fragmented systems. This is because soil macrofauna affect decomposition and also respond to the changes in abiotic conditions across habitat gradients. This study investigates edge effects on the macrofauna contributions to decomposition across forest-field ecotones. We used bait lamina assay to quantify aboveground and belowground feeding activity of soil macrofauna in Autumn 2025 in three deciduous forest-old field ecotones and one coniferous forest-old field ecotone, in Southwestern Ontario, Canada. Vegetation diversity and composition, LAI and soil characteristics (i.e., soil organic matter, pH, temperature and moisture) were measured at each plot along the ecotone. Pitfall trap data collected in Summer 2025 at the same sites were used to characterize macrofauna communities. We used generalized linear mixed effects models to estimate the effect of distance to edge, site, and depth into the soil on bait lamina consumption and soil macrofauna, with transect nested within site as random effects. Consumption activity increased with distance into the forest from the field, with the edge representing an intermediate; and, decreased with increasing depth into the soil. In contrast, soil macrofauna abundance, especially isopods, decrease with distance into the forest from the field. These trends varied significantly across sites, so that consumption activity and abundance sometimes remained constant across the ecotone (i.e., site x distance interaction). The results demonstrate that macrofaunal contributions to bait consumption varied along the ecotone, shaped by interacting environmental gradients and shifts in community composition unique to each site.

Urease is a nickel-dependent enzyme that plays an important role in urea hydrolysis and in a process named as microbial-induced calcium carbonate precipitation (MICP), which is widely used in sustainable environmental biotechnology. Despite its ecological importance, urease powers Biogrout (biocementation), a promising green technology for soil stabilization and infrastructure repair. Yet, the relationship between nickel availability, enzyme activation, and bacterial fitness remains poorly understood. In this study, we reveal a striking dual effect of nickel on Sporosarcina pasteurii: while high Ni{superscript 2} concentrations strongly inhibit growth (IC {approx} 637.7 {micro}M), they simultaneously boost specific urease activity up to six-fold. This uncoupling between biomass and enzymatic efficiency highlights a previously overlooked adaptive strategy under metal stress. Using structural bioinformatics and molecular docking, we show that Ure1--the catalytic subunit--exhibits the strongest nickel affinity (-4.3 kcal{middle dot}mol-{superscript 1}), supported by highly conserved active-site residues, whereas accessory proteins UreE and UreG display moderate and weak binding, consistent with their roles in metal delivery and GTP-dependent maturation. In addition, microscopic observations confirmed that calcium carbonate precipitation was most pronounced at intermediate nickel concentrations (approximately 400-1000 {micro}M), whereas higher concentrations ([≥]1000-1300 {micro}M) led to reduced mineral formation due to loss viable cells. Taken together, these results indicates that nickel availability controls both urease activation and bacterial fitness, and that an optimal balance is required to maximize biomenerilization efficiency in environmental applications, particularly in biocementation technology. ImportanceUrease-driven biomineralization is widely used in sustainable technologies such as soil stabilization and self-healing concrete. However, optimizing these systems requires a clear understanding of how environmental factors influence enzyme performance. This study shows that nickel, an essential cofactor for urease, plays a dual role by enhancing enzymatic activity while inhibiting bacterial growth at high concentrations. By integrating experimental data with computational analysis, we demonstrate that efficient biomineralization depends on maintaining nickel within an optimal range that balances enzyme activation and microbial viability. These findings provide practical guidance for improving biocementation processes and highlight nickel as a key regulator of urease-based environmental biotechnology applications.

Rice, a staple crop of nearly half of the world population, is grown predominantly in flooded paddies which are one of the largest contributors to methane emissions. An effective approach is to minimise the anaerobic flooded conditions that favour the growth of methanogenic archaea. Empirical measurements showed that controlled irrigation regime reduces methane emissions by 70% to 90%. The soil microbiomes of both flood and drip irrigated soil were characterised using whole-genome shotgun metagenomics. Controlled irrigation was shown to suppress methanogens and lower methane emissions. While emissions are correlated with mcrA gene abundance, empty flooded fields exhibited relatively high mcrA levels above baseline despite undetectable methane emissions. Rice cultivar genotype had no significant effect on the soil microbiomes. Co-occurrence network analysis indicates that soil microbial communities stratify according to their oxygen preferences along a gradient. Methanogens were increased in flooded paddies, and methane production attributed to the microorganisms involved in the anaerobic decay of organic matter. Controlled irrigation altered the microbiome by raising the soil redox potential by enhancing aeration and promoting ammonia oxidation and nitrification pathways. IMPORTANCEThe temporal dynamics of microbial communities in drip-irrigated rice fields remain poorly characterized to-date. Empirical measurements demonstrate that controlled drip irrigation effectively suppresses methanogens and lowers methane emissions by up to 90%. Statistical analysis further revealed a moderate correlation between methane emissions and the mcrA gene with R = 0.6 and p-value = 2.9e-05. The correlation plot showed that the outliers corresponded to samples from empty flooded fields, where high mcrA gene abundance was observed despite low methane emissions. Methane produced in the soil is likely released into the atmosphere via transport through the aerenchyma of rice plants. Controlled irrigation is shown to be climate friendly as it reduces methane emissions by improving soil aeration and increasing the soil redox potential.

The tomato leafminer, Tuta absoluta, poses a severe global agricultural threat due to its rapid leaf-mining behavior and swift development of resistance to conventional chemical pesticides. While microbial chitinases are potent biopesticides, their field efficacy is limited by environmental degradation and the short exposure window before larvae penetrate leaf tissues. This study evaluates a stimuli-responsive, controlled-release nanobiopesticide system utilizing a novel chitinase from newly characterized Serratia marcescens GBS19. A 61.1 kDa chitinase (GBS19_ChiA) was heterologously expressed in Escherichia coli and purified to a specific activity of 215.01 U/mg. The enzyme was immobilized onto starch-coated silica nanoparticles designed for target-triggered release via host alpha-amylase. Genomic profiling and R-based kinetic modeling were integrated to evaluate the efficacy of purified and immobilized forms against T. absoluta. Immobilization enhanced thermal and pH stability, with the nanocarrier maintaining 85% activity over 10 weeks. In larval bioassays, immobilization increased mortality from 21.9% to 59.4% (5000 U/mL) by day 3, reaching 62.5% by day 6. Genomic analysis identified an expansive secretome and a Type VI Secretion System (T6SS), characterizing GBS19 as a multi-pronged pathogen. Kinetic modeling established that while immobilized enzymes are effective, the 2.5-hour exposure time on T. absoluta requires the synergistic action of chitinases (ChiA/B/C) to reach the lethal desiccation threshold before larvae establish protective mines. Starch-coated silica nanoparticles significantly improve chitinase stability and delivery. However, overcoming the rapid penetration of T. absoluta necessitates a whole-cell or multi-enzyme synergistic approach to outpace larval behavioural defences.

Profiling microbiomes is an important way to understand the function and composition of communities in the wild, but natural microbiomes are often highly complex and often unamendable to experimentation to reveal cause and effect relationships. By using a small group of cultivable strains to represent those found in the wild, synthetic communities are one solution to this problem. Here we describe the MAize Rhizosphere Synthetic Community (MARSc), a genome-enabled 31-member bacterial community representative of the diversity found on the roots of maize grown in Iowa soils. This community is built around Pseudomonas putida KT2440, a model maize rhizosphere colonist and synthetic biology chassis. We characterized microbe-microbe interactions and biofilm formation of MARSc members in a variety of environmental contexts, finding that both behaviors are broadly controlled by nutrient levels. Genomic analysis and microbiome profiling of these organisms revealed that annotated biofilm genes (such as surface attachment and exopolysaccharide production) correlated to rhizosphere colonization, but neither trait correlated to in vitro biofilm formation. In vitro interactions assay findings were surprisingly consistent with co-correlations of rhizosphere abundance amongst MARSc members on roots. Finally, we found that when applied to the roots, MARSc can increase maize growth under nitrogen-limiting conditions. Altogether, MARSc is a useful tool for identifying some of the factors influencing rhizosphere microbiome assembly and will be a strong foundation for further work in this area.

Plant growth-promoting bacteria are gaining significant attention as promising biofertilizers. However, the inconsistency between in vitro plant growth-promoting traits and actual field performance remains a challenge, driven partly by a limited understanding of in situ colonization. This study characterized the colonization patterns of Citrobacter sedlakii CESi7, a novel plant growth-promoting bacterium, isolated from oil palm waste compost, during Brassica rapa cultivation. The in situ behavior of CESi7 was observed in both sterilized medium and non-sterilized soil using fluorescence in situ hybridization with a strain-targeting probe. The results revealed that CESi7 can establish both epiphytic and endophytic populations that transiently colonize roots. In a sterilized medium, CESi7 was widely distributed throughout the root tissues. Conversely, in non-sterilized soil, the bacterium formed dense aggregates specifically at the root tips. This study provides direct microscopic evidence of the colonization strategy of CESi7, offering crucial insights for its development as an effective biofertilizer.

Metabolites from Plantago lanceolata (plantain) biomass have been linked with biological nitrification inhibition (BNI) in soil. After grazing, leaf metabolite chemistry is altered via digestion, and a suite of secondary metabolites are then delivered onto soil via dung and urine. The purpose of this study was to establish if urine from sheep grazed on plantain had BNI activity when added to pasture soil, and to identify the metabolite profile(s) that most likely contribute to the BNI effects observed. Groups of sheep (n=5) were grazed on one of nine different plantain cultivars in autumn and spring with analysis of leaf material, urine, soil incubation and BNI bioassay data used to identify potential metabolite candidates implicated with BNI. The urinary nitrogen and metabolite composition of sheep fed plantain varied significantly between cultivars and season. After 28 days of incubation, all soil microcosms treated with plantain-derived urine had up to 35% less nitrate than comparative ryegrass urine controls in both seasons, except one in autumn. The key phytochemistry associated with lower soil nitrate concentrations was phenylethanoid and iridoid glycosides resulting in a higher output of glucuronidated, methylated and sulfated secondary metabolites in the urine. Among 19 secondary metabolites identified in the urine, hydroxytyrosol-related metabolites as well as catechol glucuronide, 2-methoxyphenyl sulfate and guaiacol-{beta}-D-glucuronide appear to be the most likely target compounds with respect to the BNI effects observed. Variation in metabolites from different plantain cultivars affected the ratio of metabolite derivatives in urine, which ultimately affected soil nitrification rates. Cultivar phytochemistry is therefore an important consideration with respect to BNI under urine patches. HighlightsO_LISheep grazing different plantain cultivars had different urine compositions C_LIO_LIUrines elicited biological nitrification inhibition (BNI) in soil and in vitro C_LIO_LIDifferent BNI response was related to differential expression of urine metabolites C_LIO_LIKey urine metabolites associated with BNI are derived from glycosidic compounds C_LI

Herbicide resistance threatens effective weed control in modern agriculture, particularly in grass weeds such as Alopecurus myosuroides and Lolium multiflorum. Cinmethylin is a pre-emergence herbicide with a novel mode of action that inhibits plastidial fatty acid thioesterases (FATs), enzymes essential for fatty acid biosynthesis. Although no cases of field resistance to cinmethylin have been reported, its resistance risk has not been fully assessed. In this study, we biochemically characterized defined amino acid substitutions in FAT A and FAT B to evaluate their effects on cinmethylin inhibition profile. Some substitutions in FAT A reduced inhibition in vitro, with mutations at residue R171 causing the largest shifts in sensitivity. However, these highly resistant variants required multiple specific nucleotide polymorphisms and are therefore predicted to be unlikely to arise in weed populations. In FAT B, sensitivity shifts were generally moderate. Importantly, most substitutions that reduced cinmethylin sensitivity also impaired enzymatic activity, suggesting limited viability in planta. Overall, these results indicate that while theoretical target-site resistance mechanisms exist, the practical risk of rapid resistance evolution to cinmethylin is low, supporting its value for integrated grass weed management

Root-associated microorganisms play a pivotal role in helping plants adapt to drought stress. However, the underlying mechanisms of the rhizospheric microbiome under limiting soil moisture remain largely unresolved. Integrating total and active microbiome analyses enables a more accurate interpretation of microbial responses to climate change-associated water stress. We assessed the effect of reduced rainfall on two wheat genotypes, drought-tolerant (DT) and drought-sensitive (DS), using rainout shelters that allowed 100%, 75%, 50%, and 25% of natural precipitation to reach the crop. At the peak of the growing season, rhizosphere samples were collected for metagenomic (MG) and metatranscriptome (MT) sequencing. In parallel, rhizosphere volatile organic compounds (VOCs) were collected and analysed. Differential expression analysis of metatranscriptomic data using metagenomic abundance as a cofactor was performed by comparing all treatments to the 100% precipitation control. Our results demonstrate that particularly oxidative stress-related transcripts intensify in DS as rainfall decreases. Transcriptomic shifts primarily involved upregulation of transcripts associated with antioxidant (catalase, superoxide dismutase), heat shock proteins (Hsp10, Hsp60, DnaK/DnaJ, GroEL, GroES), as well as microbial functions related to osmoregulation, proline and glycine betaine (PutA, PutP, OpuBB), and plant growth-promoting traits such as auxin production, phosphate solubilization. Moreover, volatile organic compound (VOC) emissions differed significantly between the control and drought treatments, with higher emissions, particularly acetates, in the DS genotype than in the DT genotype. Overall, pronounced drought-induced shifts in active microbial functions and VOC emissions indicate high sensitivity and functional plasticity of the active microbiome, whereas the total microbiome remains robust under medium drought.

Flavins are organic redox cofactors central to metabolism and uniquely capable of acting as extracellular electron shuttles. For life to have emerged, it must have disengaged itself from its stationary geochemical environment, a step requiring mobile redox-active components. The role of flavins at life's origin has been debated for decades, centered on their capacity for both one- and two-electron chemistry, distinguishing them from nicotinamides and iron-sulfur clusters. Here we chart the abiotic reduction of flavin mononucleotide (FMN), flavin adenine dinucleotide (FAD), and riboflavin under hydrothermal conditions (40 {degrees}C, 1 bar N2 or 5 bar H2, pH 6, 8, and 10) by nickel (Ni) and iron (Fe). Flavins show greater environmental versatility than hydride carriers such as NAD and can harvest electrons from metals that would otherwise reduce water's protons to H2. Reduction is favoured under acidic conditions, while increasing molecular charge at higher pH impedes electron transfer. Ni acts as a hydrogenation catalyst, reducing deprotonated flavins via hydride transfer, suggesting mineral composition could have influenced geochemical selection of early electron carriers. Reduced FMNH2 and FADH2 were tested as electron shuttles toward Fe3+-containing minerals, revealing that FMNH2 enables faster mineral dissolution than FADH2. We further demonstrate complete redox cycling of FMN through Ni-assisted H2 reduction and subsequent oxidation by magnetite (Fe3O4) under inert atmosphere, releasing Fe2+. This study highlights the versatility, stability and redox chemical capabilities of flavins in prebiotic context.

Mangrove forests are major blue carbon ecosystems but are often characterized by low surface methane (CH4) emissions. Such low emissions, however, do not necessarily indicate weak methanogenesis, because CH4 production may be offset by internal CH4 consumption before reaching the atmosphere. Although previous community, genomic, and transcriptomic studies have implicated methylotrophic methanogenesis in mangrove sediments, direct taxon-resolved evidence linking methylated carbon assimilation to potentially active methanogens remains limited. Here, we combined methanogenic activity assays, DNA stable isotope probing (DNA-SIP), mcrA and 16S rRNA gene analyses, and phylogenetic comparisons to identify potentially active methanogens across saline-influenced mangrove soils. The results showed that CH4 production potentials were consistently dominated by methylotrophic pathways (1.86-2.78 g CH4 g-1 soil hr-1) across all sites. DNA-SIP, together with consistent community patterns in fresh soils, indicated the potential activity of methylotrophic and mixotrophic methanogens under saline conditions. Methanolobus-affiliated methanogens were associated with salinity, Na+, Cl-, and NH4+, whereas Methanosarcina and unclassified Methanosarcinaceae were linked to soil soluble organic carbon availability and water content, indicating niche differentiation among active methanogenic groups. Phylogenetic analyses incorporating reference sequences from diverse environments further showed that potentially active mangrove methanogens were dominated by saline-associated lineages. Together with our previous methanotrophic evidence from the same sites, these findings suggest that low CH4 emissions from mangrove blue carbon ecosystems can mask substantial internal CH4 cycling sustained by active methanogenesis and CH4 consumption.